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The prediction of cross sections for nuclei far off stability is crucial in the field of nu- 
clear astrophysics. For spherical nuclei close to the dripline the statistical model (Hauser- 
Feshbach) approach is not applicable and direct contributions may dominate the cross 
sections. For neutron-rich, even-even Sn targets, we compare the resulting neutron cap- 
ture cross sections when consistently taking the input for the direct capture calculations 
from three different microscopic models. The results underline the sensitivity of cross 
sections calculated in the direct model to nuclear structure models which can lead to high 
uncertainties when lacking experimental information. 



1. INTRODUCTION 

Explosive nuclear burning in astrophysical environments produces unstable nuclei which 
again can be targets for subsequent reactions. Most of these nuclei are not accessible in 
terrestrial labs or not fully explored by experiments, yet. For the majority of unsta- 
ble nuclei the statistical model (Hauser-Feshbach) can be used to determine the cross 
sections. However, for nuclei close to the dripline the level density becomes too low to 
apply the statistical model |]J and contributions of a direct interaction mechanism (DI) 
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Figure 1. Direct-capture cross sections at 30 keV for different Sn isotopes. Levels and masses are 
calculated with models by Sharma et al. || (triangles), Moller et al. ||] (dots), and Dobaczewski 
et al. pf (squares). The lines are drawn to guide the eye. 
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Figure 2. Dependence of level energies on mass number for the even-odd isotopes 125 135 Sn 
in the RMFT model [§ (right) and for the isotopes 125 " 145 Sn in the HFB model @ (left). 
Shown are the l/2~ state (open circles), the 3/2" state (triangles) and the calculated neutron 
separation energy (full circles). The lines are drawn to guide the eye. Note the different range 
in mass numbers in the two plots. 
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may dominate the cross sections. The DI requires the detailed knowledge of energy lev- 
els (excitation energies, spins, parities), contrary to the statistical model which averages 
over resonances. Lacking experimental data, this information has to be extracted from 
microscopic nuclear-structure models. 

We compare the results for direct neutron capture (calculated in the optical model [|]) 
on the even-even isotopes 124_145 Sn with energy levels, masses, and nuclear density dis- 
tributions taken from different nuclear-structure models. The utilized structure models 
were a Hartree-Fock-Bogoliubov model (HFB) with SkP force ||, a Relativistic Mean 
Field Theory (RMFT) with the parameter set NLSH [[| and a Shell Model based on 
folded- Yukawa wave functions (FYSM) [|J. 

A similar study has already been performed for neutron-rich Pb isotopes 0]. 



2. METHOD 

The cross sections were calculated in the optical model for direct capture ||, utilizing 
optical potentials derived by the folding procedure [0]. In the folding approach the nuclear 
target density is folded with an energy- and density-dependent effective nucleon-nucleon 
interaction in order to obtain the potentials for the bound and scattering states. Only 
one open parameter A remains which accounts for the effects of antisymmetrization and is 
close to unity. The densities required for the determination of the folding potentials were 
consistently calculated from the wave functions of the respective nuclear-structure model. 
For the bound states the strength parameter A was fixed by the condition to reproduce 
the given binding energy of the captured neutron. The value of A for the scattering 
potential was adjusted to yield the same value of 425 MeV fm 3 || for the volume integral 
as determined from the experimental scattering data on stable Sn isotopes . 



In order to be able to directly compare the different models, all nuclei were assumed to 
be spherical and the spectroscopic factors were set to 1. 



3. RESULTS AND DISCUSSION 

The results of the calculations for projectiles at E c , m , = 30 keV are summarized in Fig. [I]. 
For each model we calculated the capture cross section only up to the r-process path. The 
most extreme location of the path (farthest away from the line of stability) is determined 



by neutron separation energies E n ~ 2 MeV [11]. Depending on the microscopic model, 
the path will then be located at higher or lower mass numbers A. In the case of RMFT and 
FYSM it will go through A « 132-134, for HFB the path will be shifted to considerably 
higher mass numbers A fa 142-144. (The neutron dripline is also shifted to higher masses 
in the latter model.) 

Similar effects as seen in the behavior of the Pb cross sections [[J can also be found for 
the Sn cross sections. The cross section can vary by order of magnitudes when going from 
one isotope to the next and also differ vastly between the different microscopic models. 

As the capture to low-spin states (J=l/2, 3/2) accounts for the largest contributions to 
the cross section, the results are very sensitive to the presence of bound low-spin states. 
Since the microscopic models not only yield different masses (i.e. neutron separation 
energies) but also exhibit different behaviors of the level energies with changing mass, 
"jumps" and "gaps" can be seen with some models (RMFT, FYSM), whereas others 
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(HFB) result in a smoother behavior of the capture cross sections in an isotopic chain. 
This is illustrated in Fig. |2], which shows the neutron separation energy and the excitation 
energy of the l/2~ and 3/2~ states in RMFT and HFB. As long as both states are unbound 
in the RMFT, the cross sections remain low and only jump to higher values when those 
states become bound at the shell closure. As at least the 3/2~ level is always bound in 
HFB, the cross sections show a smoother behavior. The variation in the FYSM cross 
sections can be explained in a similar way. 

4. CONCLUSION 

With this work we have underlined that the calculation of purely theoretical direct 
capture cross sections far from stability still contains a large error, even when using most 
recent nuclear- structure models. In the previously discussed case of Pb isotopes ||, the 
r-process path contains nuclei in or at the border of a region expected to be deformed, 
leading to higher level densities and thus favoring the compound nucleus mechanism. 
This is not true for neutron-rich isotopes in the Sn region, especially around the neu- 
tron magic number N = 82 where the level density becomes too low for the statistical 
model. Therefore the neutron capture cross sections have to be calculated using input 
from nuclear-structure models and will be subject to the quoted uncertainties, even when 
the different models yield similar values for other nuclear properties, such as masses. 

Similar problems may be encountered on the proton-rich side when predicting proton 
capture cross sections close to the proton dripline. 
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